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FIG. 1: (Color online) Top: Statistical distributions of 
avalanche sizes S for various injection rates (clip level C=0). 
The main plot shows Ps{S){S)"/'-'^-"'> vs S/ {Sy/'-^-°'\ and 
the inset PsiS)S" vs S/ {S)^/^'^-"\ both computed for the 
value a — 1.00 that provides the best collapses. Bottom: 
Corresponding collapses of the pdf 's of the avalanche dura- 
tions T, for the value r = 1.25. 



Planet, Santucci, and Ortm Reply: 

In Ref. [l[ we reported that both the size and duration 
of the global avalanches observed during a forced imbibi- 
tion process follow power law distributions with cut-offs, 
Px{x) = axX~"^'^Qx{x/xc)- While the exponent of the 
power law appears robust within the quoted error bars, 
the cut-off of the pdf 's depends on experimental control 
parameters such as the injection rates v. 

When adimensionalising the variables, u = x/{x), we 
observed a coUapse of the various pdf's Pu{u). Thus, 
the exponent nix should be found equal to one, as ex- 
plained in 0. Indeed, we observed clearly an average 
power law exponent a = 1.00 ± 0.06 for the avalanche 
size pdf's. However, for the avalanche duration we re- 
ported a slightly larger value within the large dispersion. 
This can be attributed to the poorer statistics for the 
avalanche duration, as we explained in [ij, affecting the 
quality of the collapse and/or the fitting procedure. This 
illustrates the difficulty of extracting accurate values of 
the power law exponents with experimental data with a 
limited statistics. 

In order to extract reliable exponents, the right pro- 
cedure -as described in [H~ is to find the power law 
exponent m* that provides the best collapse of F = 
P:r(a;)(x)"'/(2-'") as a function of X = a;/(a;)i/(2-'") . In 
order to quantify the quality of this collapse, first, we 
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FIG. 2: (Color online) Isolines of the joint distributions u vs 
w' , with a fit to the crest giving the value 7 = 1.33 ± 0.12. 



have fitted Yfu = AX"" cxp(-X/B) to the set of our 
experimental data, varying systematically the exponent 
TO and leaving A and B as free parameters. Then, we 
have computed the error e = X],Jlog(yi) — log(F/it)]^ as 
a function of m. For various velocities (clip level C = 0), 
for the size distribution, the minimum value of e (thus, 
the best collapse) is obtained for a ~ 1.00 ± 0.15, while 
for the avalanche durations we obtain r = 1.25 ± 0.25, as 
shown in Fig. 1. The insets display the collapse Pj;(x)a::™ 
as a function of x/{xy^^^~"^\ showing that the scaling 
function is very well approximated by a decaying expo- 
nential. Then it is not difficult to show that the joint 
distribution of sizes and durations can be properly an- 
alyzed using u' = S'/(S')i/(2-a) and w' = T/(T)i/(2-r)^ 

with the values of the exponents previously found. We 
show in Fig. 2 that u' oc w'"' , where 7 = 1.33 ± 0.12. 

It is important to notice that the values obtained here 
are in agreement with the original ones obtained by a di- 
rect fitting of the pdf's of the raw and dimensionless data. 
Since the actual exponents m are close to one P{x/{x)) = 
P{x) (x) ~ Pix) (x) 2^ and x/{x) ~ a;/ (x) 2^ , and due 
to both the experimental noise and limited statistics, the 
collapses previously observed were reasonably good and 
nearly indistinguishable from the present ones. 

Finally, Ref. suggests that the problem of forced- 
flow imbibition might belong to the quenched Edwards- 
Wilkinson (QEW) or equivalently the C-DP universality 
class Q, a conjecture based on a possible similarity of 
the values of the exponents. This should be taken with 
some caution, however, because forced-flow imbibition is 
a non-local process [1] , while the QEW interfacial equa- 
tion describes a local interfacial dynamics. Moreover, the 
various values of the exponents reported here are in very 
good agreement with the ones obtained from phase-field 
simulations of a non-local interfacial process 
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